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We present a model for dihadron fragmentation functions, describing the fragmentation of a quark 
into two unpolarized hadrons. We tune the parameters of our model to the output of the PYTHIA 
event generator for two-hadron semi- inclusive production in deep inelastic scattering at HERMES. 
Once the parameters of the model are fixed, we make predictions for other unknown fragmentation 

■ functions and for a single-spin asymmetry in the azimuthal distribution of tt + -k~ pairs in semi- 

inclusive deep inelastic scattering on a transversely polarized target at HERMES and COMPASS. 

, Such asymmetry could be used to measure the quark transversity distribution function. 

<N ■ 

PACS numbers: 13.87.Fh, 11.80.Et, 13.60.Hb 

<: 
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I. INTRODUCTION 

Dihadron Fragmentation Functions (DiFF) describe the probability that a quark hadronizes into two hadrons plus 
anything else, i.e. the process q — > Hi H^ X . They can appear in any process where hadronization is involved, in 
particular in lepton-lepton, lepton-hadron and hadron-hadron collisions producing final-state hadrons. They carry 
00 information that is not accessible to single-hadron fragmentation functions, but on the other hand they are more 
complex to study and to measure. 

Unpolarized DiFF were introduced for the first time by Konishi, Ukawa and Veneziano pj. Their evolution equations 
have been studied in Rcfs. and more recently reanalyzed in Refs. [3- 13- E3| • All these studies focused on the 

probability of producing two hadrons with energy fractions z\ and Z2 by integrating over the invariant mass of the 
fS produced pair. However, it is fair to say that the only experimental information related to unpolarized DiFF consists of 
invariant mass spectra of hadron pairs produced in e + e~ annihilation @, 0, , Semi-Inclusive Deep-Inelastic Scattering 
qj (SIDIS) 0, 0, E| and proton-proton collisions 0, Q, IT^L Recently it has been suggested to use DiFF as tools 
e-j ! to investigate the in- medium effects in heavy- ion collisions Iri ITU lla. Il7j . To address this and other issues, it is 
r* ■ necessary to improve our knowledge of unpolarized DiFF in vacuum. 

I DiFF can be used also for spin studies. In particular, they can act as analyzers of the spin of the fragmenting 
quark [H [H III HI |2| and the y 

can be used to study vector meson polarization [2^, 0, HH El • The definition and 
properties of all possible DiFF for two unpolarized detected hadrons have been presented in Ref. |2jj up to leading 
twist, and in Ref. [28] up to subleading twist integrated over the transverse component of the center-of-mass (cm) 
momentum of the hadron pair. Despite the wealth of observables related to polarized DiFF, experimental information 
is limited HI 13, H3. 

At present, the most important application of polarized DiFF appears to be the measurement of the quark transver- 
sity distribution in the nucleon. This function, hi, represents the probabilistic distribution of transversely polarized 
partons inside transversely polarized hadrons, and is a missing cornerstone to complete the knowledge of the leading- 
order (spin) structure of the nucleon (for a review see Ref. 32j). Being a chiral-odd function, hi needs to be combined 
with another chiral-odd soft function. The simplest possibility is to consider double-spin asymmetries in polarized 
Drell-Yan processes [33l |. This option is under investigation at BNL using high-energy polarized proton-proton colli- 
sions [3 HH an d could be studied also at GSI using polarized proton-antiproton collisions j3|| |33, 138L |3!j • 

Another possibility is to measure Single-Spin Asymmetries (SSA ) in the SIDIS production of a pion on transversely 
polarized targets. Recent data have been released using proton 0,E1 and deuteron (4^ targets. Their interpretation 
advocates the so-called Collins effect by which a leading- twist contribution to the cross section appears where 
hi is convolved with the Collins function Hi, a fragmentation functions that describes the decay probability of a 
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transversely polarized quark into a single pion. However, extracting h\ from SSA data requires the cross section to 
depend explicitly upon the transverse momentum of the detected pion with respect to the photon axis |44| . This fact 
brings in several complications, including the possible overlap of the Collins effect with other competing mechanisms 
and more complicated factorization proofs and evolution equations |4^| . 

Semi- inclusive production of two hadrons |l9t |21| offers an alternative way to access transversity, where the chiral- 
odd partner of transversity is represented by the DiFF 47] , which relates the transverse spin of the quark to the 
azimuthal orientation of the two-hadron plane. This function is at present unknown. Very recently, the HERMES 
collaboration has reported measurements of the asymmetry containing the product hiH^ |48|. The COMPASS 
collaboration has also presented analogous preliminary results |49|. In the meanwhile, the BELLE collaboration is 
planning to measure the fragmentation functions in the near future [Hoi l5l| . 

In this context, it seems of great importance to devise a way to model DiFF. From the theoretical side, this can 
help understanding what are the essential building blocks and mechanisms involved in dihadron fragmentation. It can 
also provide a guidance for fits to data and further phenomenological studies. From the experimental side, a model 
could be useful to study the effects of cuts and acceptance, to estimate the size of observables in different processes 
and kinematical regimes. Our work is not the first one in this direction plL ItH fs3| . The model presented here is 
close to the one discussed in Ref. |47j . However, for the first time we are able to fix the parameters by comparing our 
unpolarized DiFF D x with the output of the PYTHIA event generator j5^| tuned for HERMES |H3|- Then, without 
introducing extra parameters, we make predictions for the polarized DiFF and the related SSA involving the 
transversity distribution h\. 

The paper is organized as follows. In Sec. [HI we review the basic formalism of DiFF and of SIDIS cross section 
for two-hadron production. In Sec. 11111 we describe our model for the fragmentation of a quark into two unpolarized 
hadrons and give analytic results for DiFF calculated in this model. In Sec. IIVI we fix the parameters of the model 
by comparing it to the output of the PYTHIA event generator tuned for HERMES kinematics. In Sec. El we show 
numerical predictions for the DiFF and for the above-mentioned SSA in the kinematics explored by the HERMES [Isf 
and COMPASS collaborations 0|. Finally, in Sec. I VII we draw some conclusions. 



II. BASICS OF DIHADRON FRAGMENTATION FUNCTIONS 



Dihadron Fragmentation Functions are involved in the description of the fragmentation process q — > 7r + 7r~ X. The 
quark has momentum k. The two pions have masses m, = 0.140 GeV, momenta P\ and P2, respectively, and invariant 
mass Mh (considered to be much smaller than the hard scale of the process, e.g., the virtuality of the photon, Q, in 
SIDIS). We introduce the vectors Ph = Pi + P2 and R — (Pi — P-i))^. We describe a 4-vector a as [a - , a + , a x , a v ], i.e. 
in terms of its light-cone components = (a ± a?)/^J~2 and its transverse spatial components. We introduce the 
light-cone fraction z = P^ jk~ and the polar angle 9, being the angle between the direction of Pi in the pair's center 
of mass and the direction of Ph in the lab frame 55] , so that the relevant momenta can be written as 



where 
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1 Note that there is a misprint in the expressions for \R\ in Eq. (27) of Ref. 55] and in Eq. (23) of Ref. |28| . 
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and 4>r is denned later in Eq. (|15(l (see also Fig.^). It is useful to compute the scalar products 

Ph ■ R = o, (5) 

R.k=(^.-z e + M )\S\^e-ii r -Rr. P) 



2z 2M h 

Fragmentation functions are extracted from the correlation function [F 



where [M 15 



A*(z,co S e,Ml<l> R ) = ^- [ d 2 k T dk+A«(k;P h ,R) (8) 
IQMhJ k-=p-/z 

I 

Ai(k,P h ,R) tJ =J2 /^^(O^^^ (9) 

Since we are going to perform the integration over the transverse momentum kr, the Wilson lines hi can be reduced 
to unity using a light-cone gauge. 

The only fragmentation functions surviving after /^-integration are |27l l55| 



D?(z,cos0,M£) = 4 7 rTr[A 9 (z,cos6',M^(/) fl ,)7~], (10) 
e l^H H ^(z,cos9,M 2 ) =4^Tr[A*(z,cos0,A4 2 ,^)*^ 75 ]. (11) 

These functions can be expanded in the relative partial waves of the pion pair system. Truncating the expansion at 
the p-wave level we obtain [55| 

D\[z, cosd, Ml) « Dl 00 (z, M 2 h ) + D\ ol (z, M 2 ) cosO + D\ u (z, M 2 ) I (3 cos 2 9-1), (12) 
H? q (z, cos 9, Ml) » ff« « (z, Af 2 ) + ff« « (z, M 2 ) cos 9. (13) 

The fragmentation function D\ MO can receive contributions from both s and p waves, but not from the interference 
between the two, D\, i and H^ ot originate from the interference of s and p waves, D\ t u comes from polarized p waves, 
and H^ H originates from the interference of two p waves with different polarization. 

Our model can make predictions for the above fragmentation functions as well as for transverse-momentum- 
dependent fragmentation functions, which we do not consider in this Section. However, we will focus our at- 
tention mainly on the functions Di, 00 and H^ ot because of their relevance for transversity measurements in 

sidis HHUmiinl. 

Let's consider in fact the SIDIS process Ip — > I' ir + ir~ X, where I and I' are the momenta of the lepton before and 
after the scattering and q = I — I' is the momentum of the virtual photon. We consider the cross section differential 
in dM"l, d(j)R, dz, dx, dy, dips, where z, x, y are the usual scaling variables employed in SIDIS, and the azimuthal 
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angles are defined so that (see Fig. 

(q X f) (qxS) , (fx S)-q 

cos0 s = ^ — 4-~ =r> sin<fc?= v -, (14) 

\qxl\ \qxS\ 7 * \qxj\\qxS\' 

(qx I) (q x R T ) . (fx R T ) -q 

cos (f>R = =■ ■ — , svn<f>R = = — , (15) 

\qxl\ \qxR T \ \qxl\\qxR T \ 

where q = q/\q\ and Rt is the component of R perpendicular to Ph- 



2 The definition of the angles is consistent with the so-called Trento conventions l58l . 



FIG. 1: Angles involved in the measurement of the transverse single-spin asymmetry in deep-inelastic production of two hadrons 
in the current region. 



When the target is transversely polarized, we can define the following cross section combinations 



= = - £ || 1 "" 1 ;y /4 *»« + « «w § aa<». <»> 

where a is the fine structure constant, 7 = 2Mx/Q, and M is the mass of the target. These expressions are valid 
up to leading twist only. Subleading contributions are described in Ref. psj . In particular, they give rise to a term 
proportional to cos^n m d^uv an< ^ a term proportional to sin^g in da UT . Corrections at order as were partially 
studied in Ref. [4|, but further work is required. 
We can define the asymmetry amplitude 



Note that we avoided simplifying the prefactors because numerator and denominator are usually integrated separately 
over some of the variables. 



III. FRAGMENTATION FUNCTIONS IN A SPECTATOR MODEL 

We aim at describing the process q — > ir + ir~X at invariant mass Mh < 1.3 GeV. To have an idea of the prominent 
channels contributing to this process, we examined the output of the PYTHIA event generator j5^| tuned for HER- 
MES which well reproduces the measured events at HERMES. Further details concerning the event generator's 
output will be discussed in the next section. Fig. [2] shows the number of counted dihadron pairs in bins of Mh (200 
bins from 0.3 to 1.3 GeV). The total amount of events is 2667889. 

A few prominent channels contribute to this process: 

1. q — > 7r + 7r~ X±: fragmentation into an "incoherent" tt + t:~ pair that we will call, in the following, "background"; 

2. q — > p X2 — > 7r + 7r~AT2: fragmentation into a p resonance decaying into 7r + 7r~, responsible for a peak at Mh ~ 
770 MeV (14.81%); 

3. q — > lo X3 — ► 7T + 7r _ X3: fragmentation into a u) resonance decaying into 7r + 7r~, responsible for a small peak at 
Mh ~ 782 MeV (0.31%); 



3 The definition of the angles in Eqs. 1141151 is consistent with the so-called Trento conventions |5S|I and it is the origin of the minus sign 
in Eq. 1171 with respect to Eq. (43) of Ref. l55l (compare <f>R and <f>$ in Fig.^with the analogue ones in Fig. 2 of Ref. l55l *>. 
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FIG. 2: Semi-inclusive dihadron counts in bins of M h from the PYTHIA event generator tuned for HERMES ; ; 54]. The 
thick solid line represents the sum of all channels. The thin solid line represents the sum of channels 2,3, and 4 described in the 
text. The dashed line represents the sum of channels 5 and 6 (which are excluded in our model). The gray line is the difference 
between the total and the sum of all channels 2 to 6 and is assumed to represent channel 1. 



4. q — > uo — > 7r + 7r X 4 with X4 = 7r° X 1 ^. fragmentation into a u> resonance decaying into ir + ir ir° (ir° 
unobserved), responsible for a broad peak around Mh ~ 500 MeV (8.65%); 

5. q — ► J] X' 5 — > ir + ir~Xs with X$ = X X' b : fragmentation into a ?7(547) or 77' (958) decaying into ir + ir~X (X 
unobserved), responsible for a peak around Mh ~ 350 MeV (2.05%); 

6. q — * K a Xq — > ir + ir~ Xq: fragmentation into a K° resonance decaying into 7r + 7r~, responsible for a narrow peak 
at M h ~ 498 MeV (3.41%). 

On top of these, there could be the presence of two other channels: 

7. q — ► (7 Xj — ► ir + n~X 7 : fragmentation into the largely debated a resonance (see, e.g., Ref. jS^) decaying into 
7r + 7r~, which could be responsible for a very broad peak anywhere between 400 and 1200 MeV; 

8. q — > /o X% — > tt + it~ Xg: fragmentation into a /o resonance decaying into tt + tt~ , which should give rise to a peak 
at M h ~ 980 MeV, not evident in the output of PYTHIA. 

In our model, we considered only channels 1 to 6. All events not belonging to channels 2 to 6 were included in channel 
1, which then contains 70.77% of the total events. 

We work in the framework of a "spectator" model for the fragmentation process: for q — > 7r + 7r~A, the sum over 
all possible intermediate states X is replaced by an effective on-shell state - the spectator - whose quantum numbers 
are in this case the same as the initial quark and whose mass is one of the parameters of the model. In principle, 
different channels could produce spectators with different masses. Moreover, each channel could end up into more 
than one possible spectator [6(J . For sake of simplicity, here we consider just a single spectator for all channels. We 
shall denote its mass as M s and its momentum as P s . The choice of using the same spectator for all channels implies 
in particular that the fragmentation amplitudes of all channels can interfere with each other maximally. In reality, it 
is plausible that only a fraction of the total events ends up in the same spectator and can thus produce interference 
effects. 

Pions in channels 2 and 3 are obviously produced in relative p wave, since they come from the decay of a vector 
meson. In channel 4, each charged pion can be in a relative p wave with respect to the other one or to ir°, the net 
result being that there is a fraction of tt + it~ pairs that is produced in a relative s wave. In the following, we will 
neglect this fraction and assume that all charged pairs are produced in p wave; at present we don't have enough 
information to discriminate the two contributions. This assumption is most probably inadequate and would lead to 
an overestimate of the contribution of channel 4 to the final single spin asymmetry. 

We further assume that all pions in channel 1 are produced in s wave. It is possible that a fraction of the background 
events are also produced in p wave. However, such a fraction cannot be too big, as it would give rise to interference 
effects that would distort the shape of the p meson peak. It is actually known that such a distortion can indeed occur, 
but also that it is not big 0, |^| . We think that this point deserves further attention, but should not change the main 
features of our results. 
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FIG. 3: Diagrammatic representation of the correlation function A in the spectator model. 



We model the correlation function in the following way (see Fig. 03 

1 (# + m 



A«(k,P h ,R) 



F s * e 

(2tt) 4 fc 2 - to 2 ' 

/ -* 2 --4 
x ( F s e ^ +F p e 



fe 2 



(# + m) 



(19) 



fc 2 



2^<5((fc-P ft ) 2 



Isospin symmetry implies that the fragmentation correlator for u — > 7r + 7r X is the same as for d — > 7r + 7r X, 
rf — * 7r _ 7r + X, -0 — * 7T _ 7r + X. Therefore, the result for d and u quarks can be obtained from the result for u quark by 
simply changing the sign of R, i.e. changing 9 — ► 7r — 9 and — ► + 7r. From now on wc will drop the superscript 
indicating the quark flavor and calculate the fragmentation functions for u — *■ tt + t:~ X. The terms with vertex F s 
refer to the s-wave contribution, the terms with vertex F p to the p-wave contribution. The exponential form factors 
suppress the contributions from high quark virtualities |r32| |. Other possibilities can be considered, e.g., dipole form 



factors |4j],|52j], or sharp cutoffs 63] . 
Inserting Eq. (jT§|) into Eq. (JSJ), we get 



A(z, cos9, M%, cj) R ) 



\R\ 



128n 2 M h 2(l-z)P" 



d\k 7 



\F 



s|2 



( # + m) ( If, - If h + M s ) ( # + m) 



(fc 2 — m 2 ) 



2\2 



\FP\ 



F s * F p e 



F s F p * e 



(y + m) ( # - Ph + Me) g ( y + m) 
(fc 2 — m 2 ) 2 

-fc 2 ^ ( I + m) ( jj - Ph + Ms) g ( I + m) 
(fc 2 — m 2 ) 2 

£(# + m) M s ) ( # + m) 



(20) 



-fc- 1 



(fc 2 -m 2 ) 2 



with fc 2 fixed by the on-shell condition of the spectator [52|, i.e., 



|fc T | 2 



Ml Ml 
(1-*) * ' 



(21) 



The first and second lines of Eq. H20|l describe the pure s- and p-wave contributions and, as such, they are the only 
ones that can contribute to the functions Di oa , D\ u of Eq. (|12f> and H^ lt of Eq. . while the third and fourth 

lines describe the sp interference and they contribute to the functions Di j0 ; and H^ ot . 
For convenience, we introduce the function 



L\z,Ml) 



1 



i 



m . 



(22) 



This function has to be always positive for kinematical reasons. 
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We obtain the following result for the unpolarized fragmentation function 

z\R 



Di,oo(z,M£) 



\F s \ 2 e 



\&i:M h 
l-z Ml - (m + M s ) 2 



i + 2 MjMm + M s )3W 2,L» 
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The incomplete T function - typically appearing in model calculations with exponential form-factors 62] - is defined as 
r(0, z) = J z e~ l /tdt. The first term of the fragmentation function can be identified with the pure s-wave contribution, 
also called D{ ot} /4 in Ref. [55|. and the second one with the pure p-wave contribution, also called 3-0^ DO /4. 
We give here the results also for the other functions appearing in Eq. I|12l) : 



D hol (z,M 2 ) 
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16tt M h 



2Re(F s *F p )e 
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2 (m + M s ) ■ 
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(24) 



27 



D 1>u (z,MZ) = -Dl 00 (z,M*) + 



16 7T M h 
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(25) 



where 2/A 2 p = 1/A 2 + 1/A 2 . 

For the interference fragmentation function H^ ot we obtain 



H?Jz,M 2 ) = - 



T P^ r 2lm(F-F p )e^ ^ 
16 irM h v ' z 2 



z M 2 



L 2 



z + 2— 
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r o, 



2z£ 2 



(l-z) A 2 



(26) 



The interference function if-w, vanishes in our model, since it would be proportional to lm(F p * F p ). It would be 
necessary to have other sources of nontrivial phases in the amplitudes for p wave production, as the ones given by 
loop corrections [64]]. 

The vertices F s and F p are essential ingredients to reproduce the correct invariant-mass behavior and to produce 
the imaginary parts required for the generation of T-odd fragmentation functions. When no resonance is present, 
we assume the vertices to be real. When resonances are present, the vertex includes the resonance propagator and 
contains therefore real and imaginary parts. We assume in our calculation that the s-wave background is free of 
resonances and thus is purely real (at tree level). This is one of the main assumptions we make in the present work 
and has critical consequences on the invariant-mass behavior of the fragmentation functions D\^ i and H^ ot . As can 
be readily seen from Eqs. (|24|l and (|26|) . assuming F s to be real implies that D\ i is proportional to Re(.F p ) and Hf ot 
is proportional to lm(F p ). 
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At this point it is worthwhile to make a comparison with the results presented in the literature in the past. In 
Ref. [2lJ , the necessary phase difference between s and p waves was taken from ir it phase shifts in elastic scattering 
data [S3. No hypothesis was made on the invariant-mass behavior of the s and p amplitudes, i.e., on the fragmentation 
mechanism. The main assumption was that the interference pattern occurring in the fragmentation process, where 
the initial state is a quark and the final state is composed of many hadrons, is supposed to be the same as in tt tt 
scattering, where initial and final states are simply two pions. In particular, the prediction of Ref. plj changes sign 
close to the p mass. The effect is in fact proportional to the sum of the real part of the p resonance times the imaginary 
part of the a plus the real part of the a times the imaginary part of the p. Both real parts change sign close to the 
p mass. On the contrary, in our approach we neglect the contribution of the a and we take a purely real s-wave 
background, but we try to take into account in a collective manner all other ways in which pion pairs can be produced 
in a semi-inclusive fragmentation process. 

A different model prediction was also presented in Ref |47| . In that model, the s-wave amplitude was modeled as 
a sequence of two single-pion emissions and was purely real, while the p-wave amplitude contained only the p decay, 
therefore having a predictive power limited to invariant masses around the p resonance. The s — p phase difference 
followed from the presence of an imaginary part in the meson propagator. The parameters of the model were fixed 
using theoretical arguments, since no experimental input was available as a comparison. As already pointed out, the 
present work is similar to Ref |47l ] , but the s-wave amplitude is modeled by means of an effective real vertex and the 
p-wave amplitude contains contributions from the p and the u) mesons. Last but not least, the parameters of the 
model are fixed by fitting the output of the PYTHIA Monte Carlo generator, which is known to reproduce very well 
the unpolarizcd data. 

Our ansatz for the vertices is 

^ s = (27) 
_ (M 2 -M 2 p )- l r p M p (Ml - M%) - iT u M u 

~ (il' r2 1 ' ' ''' 1 T "' 1 r ' 

-ifl 



{Ml - Mff + T 2 Ml J " (Ml - Ml) 2 + VI Ml 

yJ\(M*, Ml, ml) Q(M U -m w - M h ) (28) 
4 tt T u Ml ^4 Ml ml + A(M2, M 2 , ml) ' 

where X(M 2 , Ml, ml) = [M% - (M h + m^) 2 ][Ml - (M h - m^) 2 ] and 6 denotes the unit step function. The couplings 
fsi fp, fu and f'u are parameters of the model. The first two terms of F p can be easily identified with the contributions 
of the p and the ui resonances decaying into two pions. The Lorentz structure of the resonance propagators is already 
taken into account in Eq. i|19fl . The masses and widths of the two resonances are taken from the PDG [6(|: M p = 0.776 
GeV, T p = 0.150 GeV, M w = 0.783 GeV, T w = 0.008 GeV. The details of the resonance propagators could be also 
extracted from phase-shift analyses, as done in Ref. [2]]]. I n this case the contribution of the p resonance would be 



tanc^ — i tan 2 S\ 
r p M p (l + tan 2 6{Y 



/p . " \ , ,0^ , (29) 



where S l m are the mr phase shifts for the specific (l,m) channel. However, using S l m from, e.g., Ref. |67| leads to no 
significant change compared to Eq. I|28|l . especially considering the coarse level of accuracy of our model. 

The third term in Eq. (gSJl comes from the decay of an u> resonance into three pions, of which the ir° goes undetected, 
i.e., channel 0] Let's introduce the momentum P u = Ph + Po, where Po is the momentum of tt° and P 2 — M 2 . 
Taking for oj the usual Lorentz structure of a vector meson resonance, the evaluation of the correlator A involves the 
contraction 

^ f »,«^ B _ 3 , R-Po(Vh+ ffo) _ »A , \Po\\ 

+ ^If) R "-~ $+ Ml -- $ { 1+ 3MlJ- (30) 

In the last step we took into account the fact that R has no timelike component and we have to integrate over Po- 
When performing the integration over Pq, we make a narrow- width approximation and we assume that the three 
pions are produced exactly at the lu mass (M3 = M u ). In this approximation, \Pq\ is fixed, the last term of Eq. I|30|) 
becomes negligible, the contribution of channel 0] to the vertex turns out to be purely imaginary and occurs only at 
invariant masses lower than M^ — m^ = 0.643 GeV. Abandoning the narrow- width approximation has the consequence 
of smearing the invariant-mass distribution and allowing pairs to appear at invariant masses higher than 0.643 GeV, 
as well as giving a real part to the third term of Eq. (|28l) . 

Note that we sum the three contributions in Eq. I|28|l at the amplitude level. This is the first instance where 
the assumption of equal spectators has a direct consequence, and deserves further comments. Channels [21 El and 01 
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can interfere if X2 = A3 = X4. In general, only a fraction of the total events will fulfill this requirement. On the 
contrary, since we have a single spectator for all channels this is always the case in our model. That's why we add 
up the amplitudes in Eq. (|28(l . As we shall see in the next section, the effect of these interferences in the unpolarized 
fragmentation function is in any case quite small, due to the small contribution of channel [3] and the small overlap 
between channels |5] and ^ However, a similar problem will show up also in the calculation of the function H^, with 
more relevant consequences, as we shall see in Sec. E] 

Finally, we felt the need to use z-dependent A-cutoffs to have an acceptable description of the data. We used the 
following ansatz: 

A StP = a StP z^"(l-zy-'K (31) 

The total number of parameters of the model is thus 12: 3 parameters for the form-factor cutoff A s , 3 parameters 
for the cutoff A p , the couplings f s , f p , f u and f w , the mass of the spectator, M s , and the mass of the fragmenting 
quark, m. However, in the following we shall always assume m = 0. 



IV. COMPARISON WITH PYTHIA AND PARAMETERS FIT 



In order to fix the parameters of the model, we compare it to the output of the PYTHIA event generator [53 
tuned for HERMES j&J]. The events are generated in Air. Exclusive channels are dropped. The standard HERMES 
semi-inclusive DIS cuts are applied, in particular Q 2 > 1 GeV 2 , 0.1 < y < 0.85, 0.023 < x < 0.4, W 2 > 4 GeV 2 and 
the momenta of the pions (|Pi,2|) are constrained to be larger than 1 GeV. 4 The counts per M^-bin are proportional 
to the cross section of Eq. (jl6|) times 2Mj (since the cross section in the former equation is differential in dM 2 ), 
integrated over y, x, 4>Ri <Asj an d further over 0.2 < z < 0.8. For the counts per z-bin, we integrated the cross section 
over 0.3 GeV < M h < 1.3 GeV. 

In Fig. 0] the number of counted dihadron pairs is presented binned in Mh (200 bins from 0.3 to 1.3 GeV) and z 
(200 bins from 0.2 to 0.8). From the total counts, we excluded the contributions from rj and K°, i.e., channels 5 and 
6 (see Fig. 0), because they are not relevant for our purposes. The lowest dark-gray histogram represents the sum of 
the p and uj contributions (channels 2, 3, and 4), assumed to describe the whole p-wave contribution. The light-gray 
histogram in the middle is the "background" contribution, i.e. channel 1, representing the s-wave contribution. The 
upper histogram is the sum of the other two and corresponds to the total counts minus channels 5 and 6. 

Instead of leaving all parameters of the model free, for sake of simplicity we assumed the fragmenting quark to be 
massless. We take the spectator mass to be proportional to the pair's invariant mass. The number of free parameters 
we used is then 11. 

The x 2 minimization was performed using MINUET. The \ 2 function was defined as the square of the difference 
between the expected number of events in the bin and the measured value, divided by the expected number (equivalent 
to assigning a statistical error equal to the square root of the number of events in the bin). The resulting % 2 /d.o.f 
is very high, about 25. However, we believe that the main characteristics of the Mh and z shapes of the unpolarized 
fragmentation functions are qualitatively well described. The result of the fit for the s and p wave is shown on top of 
the PYTHIA output in Fig. |t| 

The values of the parameters obtained by the fit are: 



a s = 2.60 ±0.05 GeV 2 , 
a p = 7.07 ±0.11 GeV 2 , 
f s = 1197 ±2 GeV" 1 , 
f w = 75.2 ±1.2, 



(3 S = -0.751 ± 0.008, 
p p = -0.038 ± 0.003, 
f p = 93.5 ±1.6, 
M s = (2.97 ± 0.04) M h . 



7s = -0.193 ±0.004, 
j p = -0.085 ±0.004, 
f u = 0.63 ±0.03, 



(32) 
(33) 
(34) 
(35) 



The coupling constants are fixed modulo an overall normalization factor which depends on the luminosity and is 
irrelevant for asymmetry calculations. The sign of the coupling constants is also not fixed, but the relative sign of f p , 
and f' u is (see below). 

In the p-wave channel, our model deviates significantly from the generated spectrum in the region around 0.6 GeV, 
substantially increasing the y 2 . This is due to the interference between channels 2 and 4, which is not included in the 
Monte Carlo generator. At the same time, in the s-channel the curve obtained from our model underestimates the 
data in the same region. Thus, the sum of the two curves is in good agreement with the total generated spectrum, to 



4 To perform the fit, we neglected the last cut. 
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FIG. 4: Semi-inclusive dihadron counts from the PYTHIA event generator IH tuned for HERMES and results of the fit 
(a) as a function of Mh, (b) as a function of z. Solid line: p-wave contribution; dashed line: s-wave contribution; dotted line: 
sum of the two. The contributions of the r\ and K° have been excluded. 



which the Monte Carlo generator is actually tuned. The agreement would be improved further if the contribution of 
the lo were extended at higher invariant masses by leaving the narrow-width approximation for the u> resonance and 
smearing the step function in Eq. (|28[) . Note that the interference is in this case constructive because the signs of the 
couplings f p and f' w have been taken equal. If the two couplings were taken opposite, then a destructive interference 
would take place and the model would underestimate the p-wave data at around 0.6 GeV. The agreement with the 
total spectrum would then be worsened. Also the f u coupling has been taken to have the same sign of f p to avoid 
destructive interference patterns. It is difficult with the present poor knowledge to make any conclusive statement 
about p-ij interference in semi-inclusive dihadron production. However, we can at least conclude that in our model 
the best agreement with the event generator is achieved when the three couplings f p , f u and have the same sign. 



V. PREDICTIONS FOR POLARIZED FRAGMENTATION FUNCTIONS AND TRANSVERSE-SPIN 

ASYMMETRY 

Using the parameters obtained from the fit we can plot the results for the fragmentation functions Di t u, H^ ot , and 
Di, i. The function Diji is a pure p-wave function. It depends on |F P | 2 , the modulus square of Eq. 1)28(1 . and has 
a behavior very similar to D\ OQ , the p-wave part of D\ 00 . In Fig. (a) we plot the ratio between D\^u and Di^ 00 , 
integrated separately over 0.2 < z < 0.8. In Fig. 0(b) we plot the same ratio but with the two functions multiplied 
by 2Mh and integrated over 0.3 GeV < M h < 1.3 GeV. In the same figures, the dotted lines represent the positivity 
bound 

~Di i00 <D hU <3Dl 00 . (36) 

The functions D\,ol and H^ ot arise from the interference of s and p waves, i.e. from the interferences of channels 1-2, 
1-3, and 1-4, proportional to the product (f s f p ), (f s f u ), (f s /^), respectively. Since the relative sign of f s and the 
p-wave couplings is not fixed by the fit, we can only predict these functions modulo a sign. For the plots, we assume 
that the p-wave couplings have a sign opposite to f s (as suggested by the sign of preliminary HERMES data 0). 

In Fig.|Hl(a) we plot the ratio between — \R\/M h H^ ot and £>i. OOI integrated separately over 0.2 < z < 0.8. In Fig.|fj] 
(b) we plot the same ratio but with the two functions multiplied by 2Mf, and integrated over 0.3 GeV < M/j < 1.3 GeV. 
In the same figures, the dotted lines represent the positivity bound |55| 

]| "i, <^,„ (£>?.„- I !>,„). (37, 

As is evident, there are two main contributions: 

• the interference between channel 1 (s-wave background) and the imaginary part of 2 (p resonance), with a shape 
peaked at the p mass, i.e. roughly proportional to the imaginary part of the p resonance in Eq. (|28|l : 
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FIG. 5: Model prediction for the ratio Di t u/Di t00 : (a) as a function of Mh, (b) as a function of z. The dotted lines represent 
the positivity bounds of Eq. 1361 . 



1 

0.8 
0.6 



ttj 



to; 



-0.4 



0.2 




0.4 0.6 0.8 1 

M h (GeV) 

(a) 



1.2 



m 

i 



1 

0.8 
0.6 
0.4 
0.2 



0.2 0.3 



0.4 



0.5 

z 



0.6 0.7 



0.8 



O 3 ) 



FIG. 6: Model prediction for the ratio (-\R\ H? ot )/(M H D lz0 o): (a) as a function of Mh, (b) as a function of z. The overall 
sign of _ffj ot cannot be predicted by the model and is chosen to have an agreement with the sign of preliminary HERMES 
measurements |48|l. The dotted lines represent the positivity bounds of Eq. 1371 . 



• the interference between channel 1 (s-wave background) and 4 (lu resonance decaying into three pions), with a 
shape peaked at Mh « 0.5 GeV, roughly proportional to the third (imaginary) term in Eq. (|28[) . 

The two contributions have comparable size and are large. At this point, we want to stress once more that our model 
assumptions imply that the above channels can interfere in a complete way, since the spectators X\, Xi and X4 are 
the same. As already argued before, it is likely that only a fraction of the X2 and X\ states interfere, and so does a 
(in general different) fraction of the X4 and X\ states. This could decrease the sizes of the two "peaks" of Fig. [6] (a) 
and accordingly the overall size of the curve in Fig. (b) . This is beyond the reach of our model in its present form, 
but could be a way to proceed when fitting data related to H^ ot . 

In Fig. |7|(a) we plot the ratio between —D\ t0 i and Z?i j00 , integrated separately over 0.2 < z < 0.8. In Fig. [7] (b) 
we plot the same ratio but with the two functions multiplied by 2Mh and integrated over 0.3 GeV < Mh < 1.3 GeV. 
The dotted line in Fig. {7\ (a) represents the positivity bound |55| (in the second plot the bound lies beyond the plot 
range) 



Di,oi< y|^u(^,oo + f (38) 

In this case, the function D\_ i receives basically only one contribution, namely from the interference between 
channel 1 (s-wave background) and the real part of 2 (p resonance). In fact, its shape has a sign change at the p mass 
and is roughly proportional to the real part of the p resonance in Eq. H28J) . Channel 3 is negligible as in the previous 
case, while channel 4 plays no role now because we assumed it to be purely imaginary. 
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FIG. 7: Model prediction for the ratio —Di $0 i/Di !OO : (a) as a function of Mh, (b) as a function of z. The overall sign of Di j0 i 
depends on the choice made on the sign of H^ ot . The dotted lines represent the positivity bounds of Eq. 1381 . 



Next we compute the asymmetry defined in Eq. |JTSJ, integrated over all variables but one. In Fig. 0(a) we plot 



and in Fig. [5] (b) we plot 



a*4**+*b) (M , ) - c fo,dz\R\/M h H< ot ( z ,MZ) 

Jo.2 dz D lj00 (z,M£) 



A*$**»'\z) = c Io. 5 dM h 2M h \R\/M h H< ot ( z ,Ml) ^ ^ 



fl 5 dM h 2M h D 1 , 00 (z,M%) 



where 



7T 



iof d y jit ** ^T£ffi [ 4 (*) - ^))+(^» - Km 



4 dy i::: d X [ 4 ^ + tmnm + ftm ' 

We neglected strange quark contributions. The fact that the d and u transversity distributions enter with an op- 
posite sign is due to the symmetry properties of the fragmentation functions. As already discussed in Sec. IIIII the 
fragmentation function H^ ot is the same for all quarks, but the sign of sm(tfi R + <j> s ) changes for d and u. 
In Fig. (c) we plot 



( ) Cdy [ 4 m*))+(f?i*)+rt (*))]' 

where 

, * / ° 2 8 /o.g dM h 2M fc l^/Mfc tf^ (z, Mf) 



Jo 2 <** Jo.5 2M " ^W*. M|) 



(43) 



The choices of the integrations boundaries for x, y, z and Mh are inspired by the HERMES cuts [48|]. We took 
(s = 56.2 GeV 2 ) 

x min = Max[0.023, Q 2 m J(y(s - M 2 ))} , (44) 

x max = Min[0.4, 1 - (W2. n - M 2 )/(y( S - M 2 ))} , (45) 

y min = Max[0.1, QLx/H* - M 2 )),(W^ in - M 2 )/((l - x)(s - M 2 ))} . (46) 

For the unpolarized parton distribution functions we take the parameterization of Ref. |7l| . For the transversity 
distribution function, we take the estimates of Refs. [35ll6^ . l69ll70| . The sign of the preliminary data indicates that the 
s-wave and p-wave couplings should have opposite signs and thus H^ ot should be negative. The asymmetry obtained 
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FIG. 8: v4^y ^h+^s) momen t defined in Eq. at HERMES kinematics for a proton target: (a) as a function of Mh, (b) as a 
function of z, (c) as a function of x. The different lines correspond to different models of the transversity distribution function: 
dotted line from Ref. dash-dotted line from Ref. jg^, dashed line from Ref. ^7$, solid line from Ref. |35|. 
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FIG. 9: A, 



in Wij + 0s) 



moment defined in Eq. 11180 at COMPASS kinematics for a proton target: (a) as a function of Mh, (b) as a 
function of z, (c) as a function of x. The different lines correspond to different models of the transversity distribution function: 
dotted line from Ref. |6^] , dash-dotted line from Ref. [6!| , dashed line from Ref. , solid line from Ref. |23 ■ 



from our model appears to overestimate the preliminary HERMES data by about a factor 3-4. This probably 
indicates that the model overestimates in general the effect of interferences. Apart from the overall normalization, 
the height of the bump around Mh ~ 0.5 GeV seems to be too big relative to the p peak, which is probably due to 
the fact that not all the it + tt~ pairs in channel 4 should be considered in p wave. However, in order to make more 
conclusive statements it is necessary to wait for HERMES final data. Obviously, it would be better to compare our 
model with an observable where H± ot can be isolated, e.g., in e + e~ annihilation at BELLE [Hlj . 

In Fig. [5] we plot the same asymmetry as before, but for the kinematics of the COMPASS experiment. We assumed 
the same cuts as before and change only the value of s. The size of the Mh- and z-dependent asymmetries is smaller 
than at HERMES. This is due to the sensitivity of COMPASS to lower values of x, where models predict transversity 
to be small, while the unpolarized distribution functions are big. Due to the same reason, there is a much larger 
difference among the models, as they differ substantially at low x. The asymmetries could be enhanced if the low- a; 
region is excluded from the integration. 

The COMPASS collaboration has also presented preliminary data of the above asymmetry for a deuteron tar- 
get 0| ■ We plot our prediction in Fig. ^^J 5 The different isospin structure of the target, combined with that of the 
fragmentation functions in our model, decreases the asymmetry. The x-dependent asymmetry is less than half of that 
for the proton target, while the Mh- and z-dependent asymmetries are about 10 times smaller than for the proton 
target. 



Note that the preliminary measurements of COMPASS correspond to -j4™^ h+ * s '. 
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FIG. 10: A s ^ 4 ' R+4 ' s) moment denned in Eq. ^3 at COMPASS kinematics for a deuteron target: (a) as a function of Mh, 
(b) as a function of z, (c) as a function of x. The different lines correspond to different models of the transversity distribution 
function: dotted line from Ref. |68f. dash-dotted line from Ref. [62] , dashed line from Ref. IH, solid line from Ref. 0. 



VI. CONCLUSIONS 



In this paper we presented a model for the process q — > ir + n~X at invariant mass < 1.3 GeV. We used a 
"spectator" model, where the sum over all possible intermediate states X is replaced by an effective on-shell state. 
Using this model we calculated the fragmentation functions that can be defined at leading twist when considering 
only relative s and p waves of the pion pair |55|. We obtained nonzero results for four out of five of them. 

We fixed the values of the parameters of the model by comparing the unpolarized fragmentation function D\, 00 
with the output of the PYTHIA event generator (5^| tuned for HERMES . The main characteristics of the M h 
and z shapes of Di j00 are qualitatively well described. 

We made predictions for the fragmentation functions Di n, Di a i, and H^ ot . The first one is a pure p-wave function, 
it is found to be positive, about 50% of the unpolarized fragmentation function and with peaks at the p mass and at 
around Mh ~ 0.5 GeV, where the to decaying into three pions gives a large contribution. 

The function D\ iQ i arises from the interference between s and p wave. Since in our model we assumed the s wave 
to be purely real, this function turns out to be proportional to the real part of the p wave and in particular displays 
a sign change at the p mass. The size of the function is small, in particular when integrated over the invariant mass, 
due to the sign change. Our model cannot predict the overall sign of the function. 

The function H^ ot also arises from the interference between s and p waves, but is proportional to the imaginary 
part of the p wave, i.e., it has peaks at the p mass and at around Mh ~ 0.5 GeV, due to the contribution of the 
to — > 37r channel. Its size is about 30% of the unpolarized fragmentation function. Our model cannot predict the 
overall sign of the function. 

The function H^ ot is of particular interest because in two-hadron-inclusive deep inelastic scattering off transversely 
polarized targets it gives rise to a single-spin asymmetry in combination with the transversity distribution function. 
Therefore, it could be used as an analyzer for this so far unknown distribution function. We estimated this single-spin 
asymmetry at HERMES kinematics using four different models for the transversity distribution function. We found 
the asymmetry to be of the order of 10% on average. The sign of the preliminary HERMES measurements suggests 
that H^ ot should be negative. The measurement indicates that the asymmetry in our model is about 3-4 times 
bigger than the data. This probably means that our model overestimates the effects of interferences. However, final 
experimental results are needed to make more reliable comparisons. 

For COMPASS kinematics, the enhanced sensitivity to the portion of phase space at very low x induces a reduction 
in the spin asymmetry with respect to HERMES, which can largely differ depending on the model for transversity. For 
the deuteron target, the particular isospin structure, combined with that of the fragmentation functions in our model, 
induces a further reduction such that the resulting asymmetry is much smaller than for the proton, in agreement with 
preliminary data of the COMPASS collaboration. 
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